hep-th/0411174 



Black Hole Non-Formation in the 
Matrix Model 



Joanna L. Karczmarek, Juan Maldacena and Andrew Strominger 

Abstract 

The leading classical low-energy effective actions for two-dimensional string theories 
have solutions describing the gravitational collapse of shells of matter into a black hole. It 
is shown that string loop corrections can be made arbitrarily small up to the horizon, but 
a' corrections cannot. The matrix model is used to show that typical collapsing shells do 
not form black holes in the full string theory. Rather, they backscatter out to infinity just 
before the horizon forms. The matrix model is also used to show that the naively expected 
particle production induced by the collapsing shell vanishes to leading order. This agrees 
with the string theory computation. From the point of view of the effective low energy 
field theory this result is surprising and involves a delicate cancellation between various 
terms. 
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1. Introduction 

Despite much recent progress, a number of interesting questions concerning two- 
dimensional string theories and their matrix model solutions remain unanswered. Among 
them is the seemingly simple question of whether or not the theories contain black holes. 
If they do, we would like to understand the quantum behavior of the black holes. If they 
do not, we would like to understand why not. 

The continuum 2d string theory has an exact classical solution given by an 
SL(2, R)k/U(l) worldsheet CFT [Q. For large k this CFT has an unambiguous inter- 
pretation as a spacetime black hole P,[|,|l|]0. Moreover, for the small values of k relevant 



1 More precisely for large k it is a factor of the CFT for the near extremal black fivebrane 
in critical string theory For small k, which is relevant to the present work, the black hole 
interpretation is less clear, as will be discussed in the concluding section. 
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to 2d string theory, an exact matrix model (significantly involving the non-singlet sector) 
which generalizes this classical solution to the full quantum theory is known [[5]. On the 
face of it these facts suggest that black holes indeed arise in 2d string theories much as 
in their higher dimensional cousins. Furthermore, the computations of show not only 
that the SL(2)/U(1) black hole arises in the non-singlet sector of the matrix model but 
also that there is a wider family of black holes than we would have suspected from the low 
energy effective action, since one can vary the temperature. 

On the other hand, so far there has been no clear indications of black hole forma- 
tion in any matrix model scattering amplitudes ||, (see || for an early discussion). 
There is also no sign of a large degeneracy of states which would be required for a mi- 
croscopic accounting of black hole entropy, or the phase transition signaling black hole 
formation. Finally, the quintessentially thermal nature of black holes seems at odds with 
the integrability of the matrix model H . 

Since all of these facts involve purely the singlet sector of the matrix model, all that is 
known so far is consistent with the black holes arising only in the non- singlet sector. Since 
the non-singlet sector is much harder to analyze in the Lorentzian context than the singlet 
sector, we would like to make sure we are not missing some black hole configurations that 
arise in the singlet sector. 

In this paper we address this problem by trying to make black holes by collapsing 
the tachyonic or axionic matter. It is shown that the low-energy effective action indeed 
has solutions in which such collapse occurs in a region where string loop corrections can 
be made arbitrarily small. This suggests that the problem of black hole formation is a 
question for classical string theory or equivalently worldsheet conformal field theory. On 
the other hand, a' corrections can never be controlled, so a full classical string analysis is 
needed. Using the matrix model, we show that a typical collapsing pulse does not form a 
black hole. Rather it bounces off of the would-be event horizon just before the black hole 
is about to form, and the matter is backscattered back out to infinity! 

We further analyze the problem of particle creation by the collapsing shell. Naively, 
one expects a rising tail of outgoing Hawking radiation even before a black hole forms. 
A matrix model computation is used to show that the leading term (at early times) of 
tachyon particle production on X + vanishes. A consistent picture emerges from a careful 
analysis of the spacetime effective action, which reveals several sources of particle produc- 
tion in addition to the Hawking radiation. For appropriate numerical values of the effective 
action coefficients, these different contributions can cancel, in agreement with the matrix 
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model computation. We also perform the exact (in a') string theory scattering amplitude 
computation and show that the amplitude vanishes here (this is the bulk (2, 2) amplitude 



studied in 10,11 



Our results are certainly consistent with the idea that there are no black holes in the 
singlet sector of the matrix model. At the same time, our investigations have uncovered 
more curious behavior indicating there is much we do not understand about the matrix 
model. As such we feel it is too early to jump to conclusions. Further discussion can be 
found at the end of this paper. 

This paper is organized as follows: In section 2 we describe the collapse of a matter 
shell in the 2 dimensional effective theory for OB string and discover that while other 
corrections are under control, the curvature on the would-be horizon is of order one and 
thus the effective description breaks down. In section 3 we describe the same pulse in 
matrix theory and discover that most of the energy does in fact escape to infinity at 
parametrically the same retarded time at which the horizon would have formed. In section 
4 we switch gears and begin to study early time particle creation in the collapsing pulse. 
We consider various sources of early time outgoing stress energy for both the tachyon and 
the axion in effective theory. In 4.4 we study the related string theory S-matrix element. 
In section 5, we study this outgoing stress energy in the matrix model at \i = using the 
free fermion formulation. Finally, in section 6, particle creation due to a very energetic 
incoming state is considered, with implications for possible thermality of the outgoing 
state. We close with discussion in section 7. 



Related independent work has recently appeared in |[12| . Previous work on particle 
creation in the matrix model includes EMU 



2. Spacetime gravitational collapse 

The leading bosonic terms of the low-energy effective action for type OB string theory 

are 

Seff = 7 L /^v /II ^(e- 2$ (l6 + i? + 4(V$) 2 -i(VT) 2 + 2T 2 )-i(VC) 2 + ). (2.1) 

In this expression a! = h, T is the NS-NS sector "tachyon" and C is the RR axion. The 
position-dependent string coupling is given by the dilaton 

9s = e* (2.2) 
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The corrections to ( |2.1| ) involve higher derivatives, more powers of g s and/or higher non- 
linear dependence on T. To all orders in perturbation theory there is a RR shift symmetry 
which prohibits the axion C from appearing without a derivative. The classical terms in 
the action all scale like A 2 under $ — > $ — In A, C — > AC, so that the classical limit is 

C-e-*- — ^oo. (2.3) 
9s 

In this section we study black hole formation using the effective action ( |2.1|) , and also 



discuss when and how the approximation (2.1) to the exact theory becomes unreliable. 



2.1. The tachyon 

The question of whether or not gravitational tachyon collapse occurs at c = 1 was 
discussed in the early days of the bosonic matrix model, which does not have the axion. A 
naive argument that there is no reliable approximation in which the tachyon field can be 
seen to form an event horizon goes as follows. In order for an event horizon to form, the 
metric and or dilaton must be of order one. The equations of motion following from ([2.1] ) 
allow this to happen only if the tachyon T is also of order one. However if the tachyon is 



of order one, its self interactions are important, and (|2.1| ) can not be trusted. Hence one 
cannot be sure that gravitational collapse of tachyons can really occur. 

This argument is a bit too fast because we should also consider the possibility of 
making T small by spreading out the pulse. However the problem of tachyon collapse is 
still difficult to analyze in part because of the complicated tachyon-dilaton interactions 



present in (|2.1| ), as well as the fact that the stress tensor is not positive definite even at 
leading order. We therefore turn to the problem of axion collapse which turns out to be 
much simpler. 

2.2. The axion 

Setting the tachyon to zero and scaling the axion as in ( |2.3| ), the action ( |2.1j ) becomes 
precisely the classical CGHS IT7| action coupling dilaton gravity to conformal matter.@ It 



is convenient to choose conformal gauge g^ v = e 2p rj^ l/ or, in light-cone coordinates, 

ds 2 = -e 2p dx+dx-. (2.4) 



2 This is not exactly a coincidence, since the CGHS theory was derived from the linear dilaton 
theory in the NS fivebrane throat. 
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We then have R = 8e 2p ds r d-p and the equations of motion become 

$ : e - 2 (*+p) [-4d + <9_$ + 4d + $<9_$ + 2d+d-p + Ae 2p ] = 0, 

p : e" 2 * [2d + d_$ - 4d + $d_$ - 4e 2p ] = 0, (2.5) 

c : a+a_c = o. 

Since we have gauge fixed g + + and g__ to zero we must also impose their equations of 
motion as constraints. These are 

e" 2 *(4d ± pd±$ - 2d± 2 $) = -if,, 

(2.6) 

Tf ± = \d±Cd±C . 
We are interested in the general infalling axion solution which is 

C(x + ,x~) =C(x + ). (2.7) 

Further simplifications are obtained using the residual gauge freedom to choose Kruskal 
gauge 

$ = p, (2.8) 

so that 

ds 2 = -e 2 *dx + dx~. (2.9) 

In this gauge the linear dilaton vacuum is e -2 * = — 4x + x~ with the Kruskal coordinates 
running over the range — oo < x~ < 0, < x + < oo. The ++ constraint equation is 

The general solution of (2J3) and (|2.10| ) is then 

/x + py + 
dy + / dz + T^ + (z + ). (2.11) 



Let us send in a pulse so that T?, is non- vanishing only between x + = 1 and x + = e L . 
(This is an interval of length L in the asymptotically inertial coordinate a + = lnx + .) The 
total mass of the pulse as computed on X~ is 



M = I dx + x + T° + . (2.12) 







Fig. 1: Penrose Diagram for the formation of a black hole in 2d string 
theory. 

The general solution ( |2.11p for such a pulse is depicted in figure 1. An event horizon, 
defined as the boundary of the region from which light rays can escape to Z + , forms at 
some value of x~ which we denote xjj. There is also an apparent horizon, defined as the 
line along which the dilaton gradient turns from spacelike to null so that d+e~ 2 ® = O.i 

3 In the derivation of ( (2.1|) from dimensional reduction, e~* originates as the radius of the 
higher dimensional spheres. Hence from this perspective the apparent horizon is the boundary of 
the region of trapped spheres JisR . 
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Classically the apparent horizon is always inside the event horizon. It crosses the point 
(e L , xjj), i.e. the intersection of the event horizon with the top of the pulse. It follows that 

= d + e- 2 ' i> (e L ,xJ I ) 
= / dx + d\e 2 *(x + , x H ) + d + e 2 ®(l,x H ) (2 13) 

= - J dx+T^ + (x+) + d+e-^il, x~ H ). 
Equivalently, using the fact that below the pulse we have the linear dilaton vacuum, 



-4xjj = j dx + T° + . (2.14) 
We note that the maximum value of the string coupling relevant to this computation 

A = ~ (2.15) 
occurs where the bottom of the pulse intersects the event horizon. 
2. 3. Corrections 

The first issue is whether or not semiclassical perturbation theory breaks down. In 
order that it not break down we need the string coupling e $ to be small everywhere 
outside the event horizon and in the non-vacuum region above the bottom of the pulse 
(x + > 1). Causality^ forbids a strong coupling region inside the horizon from affecting 
horizon formation. The maximum value of the string coupling in this region is given by 
( |2.15|) , so in order for perturbation theory to be good we need the integral in ( 2.14 ) to be 



large. This is obviously easy to arrange just by making T?, large. That is, if we throw 
in enough energy the black hole can be formed in a region of weak string coupling where 
quantum corrections to (|2.1|) can be made arbitrarily small. 

However we still need to worry about higher-derivative a' corrections to ( |2.1| ).l One 
example of such a classical term (consistent with (|2.3|)) is 



d 2 x^e 2 *(VC)\ (2.16) 



4 Of course it is not completely clear to what extent spacetime causality follows from the 
matrix model. 

5 For the tachyon self interaction terms such as e~ 2 *T 3 must also be considered. The absence 
of these for the axion facilitates the analysis. 
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However this term identically vanishes for the solution ( |2.7|) because we have taken d-C = 
0. To get something non-zero we need a term like 

J d 2 x^e 2 *{VC ■ V$) 4 . (2.17) 

This would give a correction term to the stress tensor of the form 

ATf + ~ e ^(d + C) A . (2.18) 



We need this to be small relative to the leading term (|2.14| ). This is also easy to arrange 



Choose constant d+C = f inside the pulse. Then we have from ( p. 14] ) and (|2.15| ) 



/= ^— (2-19) 
9HVe L - 1 



The correction term obeys 



^g± < 2g%(eL - l)f = 1 ^- T) . (2.20) 
So we need only take e L — 1 >> 1 in order for this to be a small correction. This shows 



that corrections of the type (|2.18| ) cannot prevent black hole formation in all cases. 

However there is another type of oi correction which can not be controlled: couplings 
to curvature. One such term is 

RiyCf. (2.21) 

Since R is of order one in string units at the horizon, such terms may have a large effect 
at the horizon. Hence we cannot be sure whether or not a black hole really forms in the 
full string theory. 

We may also consider how the background tachyon T affects the story. In the preceding 
discussion we have set it to zero. This is consistent with the leading equations of motion, 
so if there is a tachyon tadpole at some order it would be suppressed. However we could 
also consider the effects of the tachyon wall, i.e. a tachyon vev of the form 

^ (2.22) 



x + x~ 



This produces a wall at the region where T is of order one. One expects the axion to be 
reflected from the wall. This could be described for example by an interaction term 

J d 2 x^e T (VC) 2 , (2.23) 



which effectively imposes a reflecting boundary condition for C in the region where the 
tachyon is rapidly varying. The location of the tachyon wall can be characterized by g wa u, 
the value of the string coupling when T ~ 1. By taking 



we can arrange for the tachyon wall to be far behind the region where the event horizon 
forms. Hence all effects of the tachyon vev on events leading up to the formation of the 



Of course if we take the tachyon wall to be outside the region where the event horizon would 
form, the collapsing pulse should reflect back to infinity rather than make a black hole. 

In conclusion, due to uncontrollable a' corrections involving the curvature near the 
horizon, the question of whether or not black hole formation occurs cannot be answered 
within low energy field theory. However, since string loop corrections can be suppressed, 
it is a question for classical string theory or equivalently worldsheet conformal field theory. 
A construction of a worldsheet CFT describing black hole formation would certainly be of 
great interest. 

3. Matrix model picture 

The classical limit ( |2.3| ) in spacetime corresponds to the Fermi liquid approximation 
||19|| in the free fermion formulation of 2D string theory. In this section we describe the 
collapsing pulses of the previous section as finite perturbations of the Fermi surface and 
analyze their behavior. We will be working in the matrix model of the type OB string. 

3.1. Collapsing pulse as a Fermi sea perturbation 

The first step is to express the pulse in the standard linear dilaton Liouville theory 
coordinates 



g-wall 3> <?ff, 



(2.24) 




9 wall 



t ± In A, 



(3.1) 



where A is the matrix model eigenvalue governed by the Hamiltonianl 



H = \(d t \) 2 -\\ 2 . 



(3.2) 



For type (bosonic) theory this is a' 



\ («' = !)• 
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The asymptotic form of the dilaton is 



$ = -2 In A = t~ -t+ (3.3) 

and ds 2 — —dt + dt~. Their relation to the Kruskal coordinates x of the previous section 
is 

x ± = ±\ 2 e ±2t = ±\e ±2t± . (3.4) 
Thus, a constant q§tC = f becomes d+C = fe 2t+ over the region 

8 + C = fe 2t+ | In 2 < t+ < §ln2 + - . (3.5) 

Notice that the pulse height reaches a maximum value 

\d+C\ ~ 2fe L (3.6) 

which is exponentially large in L. It can be seen that this is generically the case. Actually 
in the matrix model it is awkward to describe pulses with J dC ^ 0. It is not hard to 
arrange for this integral to vanish. In these cases one finds that the pulse reaches values 
of order e L of either sign. 

According to P0,PT] a RR pulse corresponds to a left-right antisymmetric fluctuation 



of the two branches of the Fermi sea which we denote tjrr. The correspondence involves 
the leg-pole transform 

dVR R (t + ) = J dyK(y)d+C(t+ - y) (3.7) 



where the kernel is given by 



dkfw\-H'/ar(l/2 + ik/2)_ ik , 



with 

'7T\V 



z =\i) e ~ y ■ ( 3 - 9 ) 

This changes the shape of the pulse mostly on scales of order a' . In particular, as can be 
shown quite generally, the Fermi sea fluctuations corresponding to the spacetime pulses 
under discussion themselves have magnitude e L and are of either sign. As an illustration, 
a graph of the pulse 

d + C = fe 2t+ iln2 < t+ < iln(e L + l) 

T ( 3 - 10 ) 
8+C = -fe 2t iln(e L + l) < t+ < |ln2 + - 
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Fig. 2: Derivative d+C of the axion pulse (rilled) and its leg pole transform 
function of the inertial coordinate t + for L = 5. 

(which has J dC = 0), together with its leg-pole transform is given in figure 2. 

The incoming fluctuations live on the lower branch of the right hyperbola and on the 
upper branch of the left hyperbola. These are given by 



and 



which can be expanded as 



1 



dj+VRR 

2A ' 



(3.11) 



(3.12) 



(3.13) 



The pulse spills from right to left (and, correspondingly, from left to right), if at some 
point, d t +rj RR > g~\ u (d t +rj RR < -g~l a ). From ([Sip , both of these conditions are 
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satisfied somewhere as long as fe L > g~l u - With ( |2.19|) , this implies gu < e. L ^ 2 g wa u, 
consistent with ( |2.24| ), the condition that the tachyon wall is well behind the region where 
the event horizon forms in effective theory. Note that making the pulse wider (making L 
large) keeping gu fixed only makes this inequality easier to satisfy 

3.2. Outgoing radiation on2 + 

In this subsection we will see that, in the matrix model, the spilled-over peaks of 
the pulse (as well as the minima) are reflected off the strongly coupled region early and 
begin to arrive at X + at a retarded time which coincides parametrically with the time for 
semiclassical black hole formation. We shall further see that the majority of the pulse 
returns to T + in a relatively short time. We will assume that gu <C e L ^ 2 g wa ii, or that the 
fluctuations of the incoming pulse are large compared to g~^u- 

Let the generic bottom tip of the pulse be at 

dt+VRR = ^g-l u ( 3 - 14 ) 
This point moves along the hyperbola 

A 2 -p 2 = e, (3.15) 

according to 

A = v / ecosh(t + t )- (3.16) 
Using (|3.1| ) and (|3.4|), we find that in the past, the point is at 

x+ = -e~ 2t ° (3.17) 
8 

In the future, we have 

r 2 



X~ 



x - = _ e 2t _ _ ( 3 _ 18 ) 



In the past, the Kruskal coordinate of the infalling tip of the pulse is is x + ~ e , thus, 
using ( |2.19| ) and (pT6|) , we obtain that 



x ~ g H ~ x H . (3.19) 

A parallel computation can be performed for the spilling over peaks of the pulse, with the 
same conclusion. 
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Ordinarily one expects that retarded times for pulse reflection diverge for g waU ->0as 



the tachyon wall is taken to infinity, while (|3.19|) is independent of g wa u- What is happening 



is that the tips of the pulse are reflected at a time determined by their energy alone. For 



energetic pulses obeying the condition (|3.14| ) the value of g wa ii is not relevant. Moreover, 



the tip of the pulse emerges on X + where the black hole horizon is in semiclassical theory. 

According to ( |3.18| ), portions of the incoming pulse on trajectories with e.g. energy iep 
(either positive or negative) come out a time which is later by In 2 (in inertial coordinates) . 
Hence the bulk of the pulse will return to X + over a time period which is of order one. 
Ignoring for a moment the effects of the leg pole transform (see below) , the observer on I + 
will simply see a pulse that is reflected off of the surface where the black hole would have 
formed in the semiclassical theory. Hence the a' corrections act like a shield preventing 
the incoming pulse from ever getting inside its Schwarzchild radius. 

The semiclassical picture is distinctly different. Hawking radiation, consisting mainly 
of quanta with energies of order one in string units, begins to appear on X + at retarded 
times of order xjf ~ 1/<7h- This radiation continues until the black hole is fully evaporated, 
which takes a time of order the total energy E of the incoming pulse, assuming a temper- 
ature and greybody factor of order one. Hence the energy return takes a parametrically 
longer time in this picture than it does in the exact matrix model. (See for a related 
discussion.) 

Now we show that the matrix model picture is not substantially affected by the leg 
pole transform. The transform ( |3.7|) relates an outgoing collective field quantum of energy 
uj to that of the axion by the phase shift 

,„« ^ r(i/2 + W2) 

r(l/2-iw/2) v ; 

For uj order one in string units this implies a time delay of order the string time. For large 
uj we have from Stirling's formula 

Hence the leg pole factor leads to a time delay of order \yyuj for highly energetic quanta. 
Representing the pulse ( |375| ) by a coherent quantum state, one finds that a typical quantum 
has an energy which goes at most a power of the total energy E. Hence the time delays 
will go at most like InE, and the energy is returned over a time interval which is short 
compared to the interval E expected from the black hole picture. 
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4. Particle Production in the Spacetime Picture 

An infalling pulse of matter perturbs the quantum state of the outgoing matter, and 
leads to the quantum pair production of outgoing matter. There are several sources of 
this pair production. In this section we will compute the pair production in the spacetime 
picture, and in the next we will compare it to a matrix model computation. In addition to 
Hawking radiation, which exists for both the tachyon and the axion, we will also compute 
contributions due to a time dependent mass (for the tachyon) and a higher order interaction 
(for the axion and the tachyon). Additional outgoing particle flux arises from the T(dC) 2 
interaction, but we will not consider those in this section, as the resulting flux has a 
somewhat different form from the gravitational one. 

4-1. Particle creation by an infalling pulse 

The most familiar form of pair production is the Hawking radiation due to the mis- 
match of time coordinates on I + and I~ . This occurs for an infalling pulse even if a black 
hole is not formed, but is exponentially small at early retarded times. We wish to compute 
the leading exponential correction to the X + stress tensor for a massless particle, which 
may be either the axion or the tachyon, governed by the conformal gauge action 




(4.1) 



The metric conformal factor vanishes in the far past in the coordinates 




(4.2) 



and is given by 



P 



$ + t + -t~. 



(4.3) 



At a generic point, transforming from 



(|2.11|) one has 



2p = -ln(l 



U(t+)e 2t ~) ~U(t+)e 



(4.4) 



where 




(4.5) 



For a finite duration pulse U obeys 




(4.6) 
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The asymptotic stress tensor of the Hawking particles on Z + , as measured in inertial 
coordinates, is given by the Schwartzian of the transformation required to set p = on 
1 + . This is easily seen to be, again to leading order, 

T H r— = i[/(oo)e 2 *~. (4.7) 

4-2. A time dependent mass 

In this subsection we consider the general problem of particle production by a time- 
dependent mass term, 



— I d 2 u(d + SdS -m 2 (u + ,u-)S 2 ), (U 
2tt J 



where we take m 2 to vanish on X ± . Such a term is forbidden for the axion but is actually 
present, with a form given below, for the tachyon. In general this time dependent mass 
will lead to particle production. We will compute the corresponding two point function in 
the in vacuum 

A(s;t) = (in\S(s-)S(r)\in) > (4.9) 
to first order in m 2 on X + . Expanding in powers of m 2 

A = A + A! + --- (4.10) 

we find 

d 2 

A 1 (s;t) = -m 2 (s)A (s;t), (4.11) 



ds + ds 
where 

A (s;t) =ln(s~ -t~). (4.12) 

There is a similar equation obtained by the interchange of s and t. The solution of ( [4. 1 1| ) 
is 

du + du~ m 2 (u + , u~ ) A Q (u; t) + s <-> t, (4.13) 

-00 J — 00 

Defining the null integral of m 2 

/oo 
du + m 2 (u + ,u~), (4.14) 
-00 
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( |4.13|) can be written on X + (s + — > oo, t + — > oo) as 

Ai(oo, s~; oo, t~) = — / du~m 2 (u~)\n(u~ — 1~) — / du~m 2 (u~) ln(u~ — s - ). 

J — oo J —oo 

(4.15) 

We are particularly interested in the first order correction to the X + stress tensor, which 
is the s — > t limit of the second derivative of Ai 

d 2 

Ti—(t-) = 2 ^_^_ Ax(cx),3-;oo,r)| s _ =t _ = \d.m 2 {r). (4.16) 

Now we compute m 2 for the tachyon. In the linear dilaton vacuum, the wave equation 
for the rescaled tachyon 

S = e*T (4.17) 

is that of a massless field. However in a background with a nontrivial metric and dilaton 
the tachyon acquires a mass. To leading nontrivial order in a' the effective mass is 

m 2 (u + ,u-) = -4 + (V$) 2 - V 2 $. (4.18) 

Fixing conformal gauge and expanding 

$ = t~ - t + + <f>, (4.19) 

one finds to linear order thatl 

m 2 = -4(2p-d-<p + d + <p-d + d-<p). (4.20) 

Imposing the gauge condition p = (f), and using the solution ( |4.4[ ) for <fi and p, to leading 
order at early retarded times one has,i 

T__(t - ) = \d_m 2 {t-) = AU(oo)e 2t ~ . (4.21) 

Note that this has exactly the same functional dependence on the incoming stress tensor 
as©. 

In the calculation we have done, the origin of the two types of particle production 
seem quite different. One comes from a Bogolubov transformation on X + , while the other 
comes from an explicit time dependent interaction. This difference disappears in the less 
familiar null gauge 

ds 2 = -dx + dx~ + g ++ (dx + ) 2 . (4.22) 

In this gauge, unlike conformal gauge, the metric approaches unity on X + and there are 
explicit interactions between the metric and the massless field S. 



7 We are using here the conventional form of the effective action which does not contain for 
example a T 2 e~ 2 *c?-(9+ < l > term. 

8 A further correction to the mass might arise if there as an e _2 *i?r 2 term in the effective 
action. This would change the coefficient but not the functional form of the action. 
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4-3. The Quartic Interaction 

In this section we study an interaction of the form given in ( |2.16D 



^ J d 2 u(d+Cd-C + Je 2 *- 2p d+Cd+Cd-Cd-C) , (4.23) 



where J is a constant. Unlike the mass term in (|4.8Q , such a term is not forbidden for the 
axion. There is also an identical term for the rescaled tachyon field S. 

We are interested in how the incoming stress tensor T ++ affects the outgoing <9_c two 
point function. This is governed by the equation of motion 

8+d-C + 2JT++8- (e 2t ~- 2t+ d-C) = 0, (4.24) 

where we have used (fOj). Repeating the steps of the previous section, one finds the last 
term in (|4.24| ) implies an additional correction to the two-point function 



f s 2 1 

Aj(cx), s~; oo,* - ) = / du~J{u~)[ - _ _ - _ _ , 2 +] + t' 

J — oo \ U t ) \ U t ) 



(4.25) 



where 



J{vT) = -2Je 2u I du + e- 2u+ T ++ = -4Je 2u ~ U(oo). (4.26) 

J — oo 

The stress tensor correction is then 

T J __(0 = ^ F A 1 (oo,s-;oo,0| s _ =t _ =-e 2t "rj(oo)^- (4.27) 
We note that this exactly cancels the Hawking particle production ( |4.7| ) for J = 

4-4- Particle production in string theory 

We will now perform the exact string theory computation of particle production. Note 
that the process we are interested in occurs in the bulk of the two dimensional spacetime, 
far from the Liouville wall. Such processes are fairly easy to compute in string theory, 
since we can simply use the linear dilaton conformal field theory. Such amplitudes were 
computed in |I0| , [|TTf| . In order to make sense of these amplitudes one needs to consider 



oo 



wavepackets that are very localized, as explained in ||23|| . For example, we can consider 
localized gaussian wavepackets. We are interested in a process where an incoming left- 
moving RR axion creates a gravitational field which in turn creates two RR axions. So the 
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(a): (2,2) 



(b): (1,3) 



Fig. 3: (a) Diagram representing creation of Hawking radiation. The gravi- 
ton is not a physical propagating field, we are just thinking of it as an off 
shell particle, (b) Diagram used for computation of the static field created 



by the incoming pulse, as studied in ]53J 



net process is one in which we have an incoming left-moving RR axion and an out going 
left-moving RR axion and two right-moving RR "Hawking" particles (see figure 3(a)). 
Vertex operators involve a term of the form e^ - ^*. The mass shell condition isi 

(P + 2) 2 +u 2 = 4n (4.28) 

where n is an integer in the NS-NS sector and an even integer in the RR sector. The 
tachyon and RR axion have n = 0, the graviton has n = 1, etc. Massless particles have 
f3 = —2±iui, where +(— ) correspond to right (left)-moving particles. Note that left-moving 
particles are going into the strong coupling region, while right-moving particles are going 
away from the strong coupling region. Once we choose the chiralities of the particles we can 
analytically continue in 10. But we cannot analytically continue from one chirality to the 
other. Another important feature of the bulk amplitudes is that they preserve left-moving 
worldsheet fermion number (—1)^. The tachyon and right-moving RR fields are odd and 
the graviton and left-moving RR field are even. 

We are interested in a (2, 2) bulk amplitude that contains two left-moving and two 
right-moving particles. It was argued in []n[ that the (n/,n r ) amplitude (with n\ left 
movers and n r right movers) vanishes if n r > 1 and n\ > 1. Let us see how this happens 
in more detail. A direct evaluation of the (2, 2) bulk amplitude of four RR particles gives 

A XI^T R _l a\ r«/2)r(«/2)r(l+a/2) 

A.,* = S(E ft + 4)^ 2 ^4 r(1 _ t/2)r(1 _ M/2)ra _ s/2) (4.29) 



We are still in units with a' = \. 
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where we have defined 

« = = ~t[( W 3 + ^4) 2 + (ft + 04 + 2) 2 ] 

4 (4.30) 
t =same with 34 — > 31 , u = same with 34 — ► 32 

The S functions in ( |4.29| ) enforce energy-momentum "conservation". The on-shell condi- 
tions for the RR axion 

uf + (A + 2) 2 = (4.31) 

then imply 

s + t + u = l. (4.32) 

If we take particles 1 and 2 to be left-moving and particles 3 and 4 to be right-moving, 
then we have 

01,2 = -2 - iw 1>2 , 03,4 = -2 + zw 3 ,4- (4.33) 

Using the on-shell conditions (f4.31| ) one then finds 

s = 1 (4.34) 

This in turn implies that (|4.29Q is zero. 

We wish to understand the cancellations which lead to this result from the perspective 
of field theory diagrams. This can be done by relaxing the on-shell conditions ( |4.31|) , and 



defining the off shell amplitudes by analytic continuation in s, t and u in ( |4.29|) . It is 
convenient to do this in such a way that ( |4.32|) is maintained. Lets first look at the t 
channel. For s ^ 1 there are poles at t = — 2n, n > 0. These poles correspond to 
intermediate NS-NS states with (— 1)-^ = —1. The first pole corresponds to the tachyon 
and it has the simple form 

At~o • (4.35) 

This agrees with the answer we expect from the effective field theory and the coupling 
( p.23f) . Note that (|4.35|) does not vanish when we set s = 1, but all higher poles do vanish 



when we set s = 1. However, there is a similar pole in the u channel, which has the same 
form as (|4.35| ) with t — > u. Then Q4.34Q implies that u = — t so that the two diagrams 



containing tachyon exchange vanish. This also happens in the effective field theory since 
we have identical kinematic restrictions. 

A different situation arises if we look at the poles in the s channel. Notice that the 
expression (|4.29| ) has poles at s = — 1 — 2n. These poles correspond to NS-NS states with 
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(— = l. The first state is the graviton. Note that the tachyon does not get exchanged 
since the three point coupling vanishes due to the conservation of worldsheet left fermion 
numbei0. In other words, the coupling between two left-moving RR axions and a tachyon 
vanishes. The pole at s = — 1 gives a contribution of the form 

A s — i ~ — — (4.36) 

S ] L 

as in figure 3a. Note that when we set (|4.34| ) , ( |4.36| ) gives the contribution due to graviton 
exchange which can be interpreted as particle production. Since in the end the (2, 2) 
amplitude vanishes, this contribution ends up canceling after including the exchange of all 
higher massive states. One might imagine matching ( |4.29| ) to a low-energy action which 
includes only the tachyon RR axion and graviton-dilaton. In this case we would simply 
include the poles due to (|4.35| ) and (|4.36|) , then we could think about the rest as the 
quartic interaction of the previous subsection (the J term in ( |4. 23| ) ) with a coefficient 
simply given by minus ( f4.36|) with s = 1. 

There are other diagrams that could lead to particles on I + . For example one of 
the incoming (left moving) C particles can scatter into a left moving tachyon and a right 
moving C particle that goes to 1 + . This process can be explicitly subtracted off if we are 
interested in isolating the piece that comes from gravitational particle creation. In fact, we 
will see that the matrix model computation contains a piece that comes from this process. 

Note that in a rather similar diagram, involving an intermediate graviton, was 
computed. They sent in a pulse and then they send in a second pulse and they studied the 
scattering of the second pulse from the gravitational field created by the first. In this case 
the corresponding bulk amplitude is depicted in figure 3 (b). This is involves a (3, 1) bulk 
amplitude which does not vanish Jl(| . The authors of [[23| found a detailed agreement 
between the effective action computations and the matrix model results, including a non 
vanishing (1,3) amplitude. Note that the graviton which is exchanged in the case studied 
in [23] can be on shell. In other words the corresponding invariant is not fixed, but is a 



function of the external momenta. 

It seems rather surprising that the effects of static gravitational fields are seen ([ 
but time dependent ones are delicately canceled. This might be a consequence of the 
symmetry, and it might be ultimately be the reason that black holes do not form in the 
singlet sector 0. 



10 In the effective field theory the couplings coming from ( p. 23 ) vanish if the two RR axions 
have the same chirality. 
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See [24] [25] [26] [27] for further discussion on the implications of the Woo symmetry. 
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5. Particle Production in the Matrix Model 



In the previous section, we have computed particle production arising from various 
terms in the action. These results should be compared with the matrix model computation. 
The matrix model has a nontrivial S-matrix, which is further modified by the leg-pole 
factors and then interpreted as the worldsheet S-matrix of string theory. In principle, the 
various terms in the effective action which we have studied in the previous section can be 
computed from this information. An example of such computation was presented in |2ij| j, 
where the coefficients of the Einstein-dilaton term and the cubic tachyon self interaction 
where computed. 

In the matrix model, particle production at the level at which we are studying it is 
related to the 1 — > 3 S-matrix element. This should not be confused with the statement 
we made above, which was that the particle production came from a (2, 2) bulk ampli- 
tude. These two statements are consistent since one of the particles involved in the bulk 
computation then scatters from the wall and comes back out and becomes a right-moving 
particle also. The matrix model takes the full process into account. The computation we 
will eventually do is really an inclusive computation where we compute 1 — » 2 + X where 
X is anything that comes out late and the 2 represent the two right-moving particles that 
are created at early times. This represents an incoming particle causing the production of 
a particle-antiparticle pair on X + . The leading term in string coupling perturbation theory 
coming from the matrix model should reproduce the sum of the various terms computed 
in the previous section. 

In ||23|| , the computation starts out as a perturbation expansion in - or equivalently 
the value of the string coupling at the tachyon wall. At the end however, we are interested 
in the //-independent terms, which are a bit tricky to extract from a — expansion. In the 
bosonic string considered in [p3 [, one does not really have the option of simply working 
at fx = because the theory is ill-defined there. In contrast for the OB theory there is no 
problem at \i = 0, so we will compute in the relatively simple free-fermion formalism at 
\i = 0. 



5.1. Basic formulae 

In this subsection we will summarize the basic formulae needed for the computation. 
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First let us summarize the leg-pole factors. 

<wM =/(w) _1 aout(w) 
a in (<jj) =l(u)a in {to) 
T{iuj) 



bos 



r H w ) (5.1) 
Uu>\ - T(lUj/2) 



T(-iu/2) 

1{uj)c = 



Y{\ + iu/2) 



r(J-iw/2) 

where a(oj) is a creation operator for 10 > 0, overbars denote matrix model fields, T denotes 
dressed OB tachyon and C denotes OB axion. We have chosen a' = 1 for the bosonic string 
and a' = | for the superstring. The OB matrix model tachyon and RR axion are related 
to the free fermions by 

=e i(T+C) ) = jiT-C) 

_ . . (5.2) 

8T =^ e + , dC = ij e tp° + 

r ^{ e i0 + e -ic ) ^ r „ e it {e ic _ e -ic ) 

ip L (ip R ) is localized in the left (right) part of the Fermi sea, and we use ip to denote ip^ 
to avoid a cluttering of indices. Now we consider the matrix model reflection amplitudes 
for even and odd fermions: ifj e , The reflection amplitudes of fermions of energy e are 

r in {e) =ie- i/4 2-^ lf~f r out (e) (5.3) 

i bosr \ _ / r(| - ze) , 

where ^(e) is a fermion creation operator if e > 0. We have also indicated the reflection 
factor in the bosonic case (this will make sense only for appropriate Fermi levels). 

5.2. Single fermion reflection 

Let us consider an incoming fermion with a wavefunction fi n (t + x), such that fi n (t) = 
for t < 

fin(t + X ) = ^J duf in {u)e-^ t+x) (5.4) 
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Then fi n (uj) is analytic in the upper half plane. The out wavefunction is 

fout(t - x) = [ dwR{u)fMe- iu}{t - x) , R(u) = 1 { ( K ~ IUJ/ ^ (5.5) 

An important property of the fermion reflection factor is that it is analytic in the upper 
half plane. Naively this would imply that f ou t(t — x) is zero for t — x < 0. This is not 
true because R —>■ oo when uj —>■ zoo. If fact, let us consider the case where we have a 
Fermi level given by /i and we consider energies near the Fermi level. Then we just replace 
uj — > + uj in the reflection factor ( |5.5| ). For large \fj,\ this reflection factor has a leading 
phase going like R ~ e _ltjlog ' M '. This means that a low energy excitation, with energy 
uj <C will lead to 

fout(t-x) ~ fin(t-X + log\/J,\) (5.6) 

So this looks as if the reflection is happening at 2x = log \ fi\. This agrees with the intuition 
that the reflection occurs at — A 2 /2 = /i, where A ~ e x . This is a nice check of the sign in 



the argument of the gamma functions (|5.3| ) 



In other words, the early time reflection that we get in the matrix model is due to the 
growth of R in the upper half plane and not from poles in R. On the other hand the late 
time behavior of the amplitude does indeed depend on the poles in R. 

5. 3. One point functions 

We now consider an incoming coherent state 

|/) = e < / dw ^ <w)a<B(w) |0). (5.7) 

where f*(uj) = f(—uj) is the pulse profile in the matrix model and <5; n is a collective field 
oscillator.0 We want to compute 

{f\d<Kt)\f) = (W))f (5-9) 



12 These obey 



n ( 



(uj) = iu)v2 / dse lu,s rj(s), [a(u>), a(u>')] = 27tuj' 8(uj + uj'). (5-8) 



23 



Since we are interested in early time behavior on I + , we will simplify things a bit by taking 
/ not to have an early time tail and also obey J due~ u f (u) = 0. In this case the matrix 
model / does not contain the leading exponential piece at early times. 

We can compute ( |5.9|) by expressing the boson in terms of the matrix model fermions. 
For the bosonic string we obtain the expression 



(d<f>{t)) f ~i J dudee- lujt dudu'e l{Lj/2+e)u+l(uj/2 - £)u ' 



i{f(u)-f(u')) _ 1 ( 5 - 10 ) 

l(u)R-\-u/2 + e)R{u/2 + e) — 

XL Hi 

Where R{uS) are the fermion reflection factors ( |5.3|) and the symbol ~ here an hereafter 
means that factors of 2, -1 and tt are omitted. The subtraction of —1 corresponds to 
normal ordering the fermions. We are interested in the early time tail, so we want to 
consider t <C 0. In this case we may deform the ui integral to the upper half plane. Since 
the reflection factors are analytic in the upper half plane we see that (if we do not shift 
contours in e) we do not get poles from the bounce factors^!. We do get poles from the 
leg factor. For the bosonic string we would get poles at u = in, n > 0, and the leading 
contribution comes from the first one. So we evaluate the residue of the integrand at u = i. 
In the bosonic string case the bounce factors give something of the form 

R~ X R ~ |e| -> \n\ - e (5.11) 

where the second result is the one with nonzero //, \i < 0. We now set u = u c + r/2, 
u' = u c — r/2. The integrand contains a factor e ier , so the factor of e may be converted 
into id r , yielding 



e 



if(u c +r/2)-if(u c -r/2) _ 



*(H - ldr)~ -|r=0 = (W + W?) (5-12) 



r 



In conclusion, the early time piece is proportional to 

(00)/~e« [ due u (W' -\{f) 2 ) (5.13) 



where we have set u c = u. This answer is valid to all orders in perturbation theory. This 
expression is convergent if we choose / so that the first tail of / vanishes. In other words, 



13 Note that in this case it is OK to shift the contour to lo — > ioo since the growth of the bounce 
factors is canceled by the decrease of the leg factor (away from the purely imaginary axis). 
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we want / <C e~ u for u — > oo. In general, if we start with a completely localized /, then 
we get / ~ e~ nu from the leg poles. A nice derivation of this result was given in |f28[ , by 
using the fact that the early time tails are given by conserved charges in the matrix 
model. Finally note that ( |5.13|) can be interpreted as arising from a bulk interaction of 
the form je'^T 3 . 

We now compute similar expressions for the OB case. Let us start with the tachyon. 
Then we find 

(dT(t)) f = J dujdee- lujt dudu'e^ /2+e)u+l{uj/2 - e)u 'lT{oj)x 

~Rj 1 {-u/2 + e)R e (cj/2 + e)(: 4>%u')ip e (u) :) f (5-14) 
+R~ 1 (-u/2 + e)R (u/2 + e) (: : ) / 

As before, we are interested in the upper half ui plane. Only the poles of the leg factor 
contribute. The first is at uj = 2i. Then the bounce factors give 

RZ l R e = Ro'Ro = \{\ + e 2 ) = \{\ - d 2 r ) , (5.15) 

The first equality in ( |5.15| ) holds at all poles of the leg factor, namely it holds for u = i2n, 
n > and n integer. Since the two bounce factor ratios are the same we get the combination 

(: ^ e {u')^ e {u) : + : ^°{u')^°{u) :)/ = (: ^ L (u')^ L (u) : + : i> R (u')^ R (u) :) f 

e i(T («)-?(«')) cos (C(u) - C(u')) - 1 ( 5 - 16 ) 

= 1 / 

u — u 

In this and similar expressions below, T and C to denote their coherent state expectation 
values in the coherent state |/) rather than quantum operators. Hopefully the distinction 
will be clear form the context. The final result at \i = is 

(dT) ~ e 2t [ due- 2u [f + A(C') 2 f + -(T'f - -T'"] (5.17) 
J 3 3 

Let us comment briefly on this result. The term involving T 3 can arise from a bulk e - 2 *7" 4 
interaction. Notice that the difference between the bosonic string result and the superstring 
result comes from the interplay of the matrix model amplitudes and the leg factors. If we 
choose a different values of a' in the bosonic and in the superstring so that the matrix 
model amplitudes are essentially the same, then the leg pole factors are different. This 
difference in the leg-pole factors then translate into the rather different effective actions 
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in the bosonic string versus the superstring. For example the bosonic string contains a 
T 3 interaction which is responsible for ( |5.13|) , while no such interaction is present in the 
superstring. Indeed, the superstring result (|5.17| ) is related to T 4 interactions in the bulk. 



If we wanted to consider the result at non-zero /i, then we must shift e in ( |5.15|) by /i. 
Then we get terms that are at most quadratic in \i. 

Now let us consider the same computation for the RR axion 

(dC(t)) f = J du}dee- lujt dudu'e l ^l 2+e)u+l ^/ 2 -^ u 'l c {oo)* 

~R-\-u;/2 + e)R e (u;/2 + e)(:r(u')r(u) :) ( 5 - 18 ) 

+R-\-lo/2 + e)R (u/2 + e) (: */i ■) / 

As before, we are interested in the upper half u plane, and only the poles of the leg-pole 
factor contribute. The first one is at u> = i. The bounce factors give 

R~ x R e = R- X R = | (5.19) 

The first equality holds at all poles of the leg factor, namely at u = i(l + 2n). (Note that 
we get basically the same as in the bosonic string (|5.11|) ) We see that the independent 
factors then combine into 

(: 4> e (u')^°(u) : + : ^ e (u')^°(u) :)/ =(: ^ R {u')^ R {u) : - : ^ L {u')^ L {u) :) f = 

e i{T{u)-T{u')) s i n (C(u) - C(v!)) ^ 5 - 2 °) 
% 

u — v! 

Transforming (|5.19|) into a derivative we find that 



(dC) f ~ J due~ u dTdC (5.21) 

This is consistent with a bulk coupling of the form T(VC) 2 (see ( j2.23| )). 

5.4- Axion two point functions 

We now compute the connected part of the axion two point function 

{dC(t)8C(t)) f , c = 

g — i(aj 1 +LU 2 )t gi^i (ui+u[)/2+ie 1 (ui — u[ ) ^iuJ2 (u 2 +u' 2 )/2+ie 2 ("2 — u' 2 ) (5.22) 
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The ipout operators in this expression need to be replaced by in operators using the bounce 
factors as above. However, we are eventually going to be interested in in setting ui\ = i 
and u>2 = i- So in order not to get a terrible mess, we are going to evaluate just this term 
in the full answer. As we saw above taking ui\ = u>2 = i will imply that the bounce factors 
reduce to simple terms involving ie\ and ie 2 - This, in turn become derivatives with respect 
to r% = u\ — u[ and T2 = u 2 — u 2 . So we get an expression of the form 

(dC(t)dC(t)) f>c = e 2t J du c x du c 2 e-<- u * 

d ri d r2 (: $ R {u' x )^ R {ui) -^Vi^M : (5 ' 23) 

After taking derivatives, this expectation value is proportional to 

K _ 1 ^ )2 ^K)^K) 1 (5-24) 

where for simplicity we have taken an antisymmetric incoming pulse with T = 0. 
The final answer is 

(: dC(t)dC(t) :) f ~ e 2t J du x du 2 e~ Ux ~ U2 -r ^y dC( Ul )dC(u 2 ) (5.25) 

In treating the double pole we note that it arises from the contraction 

i 

%l)*{ui)%l)(u2)- : i(j*(ui)iI)(u2) :~ ~ ifj(ui)tf;*(u2)- : ^(ui)tf;* (u 2 ) : (5.26) 

ui — u 2 + ze 

Since the denominator in ( |5.25| ) came from performing these contractions we should replace 
it by l/(ui — U2 + ie) 2 . It now easy to see that there is no divergence. 

We see that ( |5.25|) does not have the interpretation of the Hawking radiation that we 



were expecting. Nevertheless we would like to identify the spacetime origin of ( |5.25| ). In 
the spacetime action we have a coupling of the form ( |2.23|) 

~ J e $ T(aC) 2 (5.27) 

This is the same coupling that was responsible for the one point function ( |5.21|) . This 
coupling can lead to a process where one incoming RR axion becomes an ingoing tachyon 
and an outgoing RR axion. By summing over the tachyon final states, then we can get 
flOp . In fact the whole answer is accounted for in this fashion. This implies that all other 
processes cancel. In fact, this is consistent with the fact that the (2, 2) bulk scattering string 
amplitude is zero. 
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5. 5. Tachyon two point function 

The situation is rather different for the two point function of the OB tachyon on X + . 
For the axion, the exponential tail was produced by the leg pole transform. The falloff rate 
had the right strength to correspond to Hawking radiation but, as we saw, the detailed 
functional form indicated an alternate source. For the tachyon, the leg pole transform falls 
of at twice the rate, so that 

(dT(t)dT(t)) f ~ e 4t . (5.28) 

Hence there is no hope of getting anything large enough to correspond to the rising tail of 
Hawking radiation. As mentioned before, this is not a contradiction because the leading 
tail of Hawking radiation can be canceled by other interactions if they have certain finely- 
tuned coefficients. In fact the bulk scattering amplitude to 2 RR fields into 2 tachyons or 



2 tachyons into two 2 tachyons vanishes [1 I 



6. Scattering amplitudes for very energetic states 

In this section we consider a single incoming tachyon in the type OA theory with very 
large energy e> 1. For simplicity we set \i = 0. We are then interested in seeing if there 
is some sign of black hole formation. Specifically we will compute the number of particles 
with energies uj that we have in the final state. This computation is very similar to those in 
0]. In particular, the analysis of is more detailed than what we are going to perform 
here. We decided to include this because is it is a rather simple and direct computation. 

The exact fermion scattering amplitude is given by0 

T(- - i-) 

™ - w4 (fU) 

An incoming boson with energy e may be written as 

\y)=a in (-e)\0)= [ dx$ n (x)ij>i n (e-x)\0) (6.2) 



In terms of out operators 

dxR*(-x)R(e-x)ijjl ut (x)^ out (e-x)\0) (6.3) 

Jo 



14 Up to phases that are constant or linear in the energy. Such terms do no matter for this 
computation. 
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We can compute the norm of fl6.3|) and we find that it is equal to 

(W|*) = e<y(0) (6.4) 

where the 6(0) is just the usual volume factor due to the fact that we are considering plane 
waves. Now we want to compute the expectation value for the number of particles with 
fixed energy u in the outgoing state (|Q|). Namely we are interested in computing 



wiVu = (^\a out (-u)a out (u)\^) . (6.5) 

Note that leg pole factors cancel out in this computation. In order to obtain N w we first 
compute that 



a out {uj)Y&) = dxF(x)ip\ ut (x - uj)ip out {e - x)|0) 

Juj (6.6) 

F(x) =R(-x)*R(e -x)- R*(-(x - u))R(e - x + u) 

where we assumed that w^e (otherwise there is an additional term) . Note «o annihilates 
the state ( |6.2| ), since ceo is the fermion number operator and the net fermion number of 
0|) is zero. 

Now we go back to ( |6.5| ) which is equal to the norm of the state (|6.6|). We find that 
) becomes (up to a 5(0)) 

dx \R*(-x)R(e -x)- R*(-(x - u))R(e - x + uj)\ 2 = 

! dx\l - R(-x)R*(6-x)R*(-(x-uj))R(e-x + u)\ 2 = (6.7) 

J U) 

I dx(2-2Re(RR*R*R)) 

where we used that \R\ = 1. Using Stirling's formula we can show that the leading phase 
at large x is 



T(- - i-' 



2 2 / „ — ixlogx 
1 i -a;' 



r(i + z§ 

We then see that for large x 

R(-x)R*{-{x -uj)) ~ e iL0losx (6.9) 
Similarly we obtain that for large e — x 

R*(e-x)R(e-x + u) ~ e -^^-^) ( 6-10 ) 
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We conclude that the product of all phases is 



x luJ (e-x)- tuJ =t luJ (l-ty tuJ , where t = ^ (6.11) 

We now substitute ( |6.11|) in the last term of the integral and change the inte- 
gration variable to t. We need to evaluate an integral of the form 

/ dtf"(l -t)-™ =T(1 -w)T(l+iu) + o(u/e) . (6.12) 

Jw/e ' SinhTTW 

Putting all this together we find that (|6.5|) has the form 



= WM^MW = 2(1 _ ™ j for u <<; e (613) 

Sinh7TCJ 

Notice that this result goes to zero for u — > as expected. 



The expectation value (|6.13|) goes to a constant for large u. If the physical process 
were black hole formation we would have expected a thermal factor with a temperature of 
order one. We might imagine that the first, constant, factor represents a reflected pulse 
that arises at the moment we form a black hole and that the second factor is proportional 
to Hawking radiation. Indeed the second factor has an exponential decay reminiscent of 
Hawking radiation. But the second factor has the wrong sign and its integral is independent 
of e@. The expectation value (|6.13| ) is saying that we have an energy of order one in each 
mode of frequency uj (for 1 C w < e). Hence there is no indication from this computation 
that a highly energetic tachyon forms a black hole. Similar observations were made in 



7. Discussion 

We close with a discussion of several issues concerning black holes in the matrix model 
in the context of our results. 

15 By writing sinl ^ VUJ ~ e~ nuj ~ for large uj in fl6.13|) we can read off the a temperature 

temperature which is the same as the temperature of the OA SL(2)/U(1) black hole, = 2ny/a'/2. 
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7.1. The definition of time 

One issue concerns the time coordinate to be used in the spacetime/matrix model 
dictionary. In spacetime computations, there are several natural retarded time coordinates 
on X + . One of them is where the conformal factor p asymptotically vanishes and the linear 
dilaton is time-independent. We refer to these as inertial coordinates, as they are associated 
to asymptotically inertial observers. The null coordinates 

ds2 = ~ dt tnertial dt inertial + 9++ ( dt tnertial ) ^ C 7 ' 1 ) 

can be seen to be inertial in this sense. 

A more common coordinate choice (employed for example in this paper and |f23| ) is 
conformal gauge 

ds 2 = -e 2p dt+ f ,dt~ f , . (7.2) 

conformal con formal \ I 

The conformal retarded time and inertial retarded times are related at early times on X + 
via 

Conformal = inertial ~ |^(oo)e 2i — « + • • • (7.3) 

These two coordinates systems define outgoing vacua on X + which are related by a non- 
trivial Bogolubov transformation. If we work in the conformal gauge, Hawking particle 
production is seen indirectly as a consequence of this Bogolubov transformation to inertial 
coordinates. In the inertial coordinates (|7.1| ), particle production arises directly from the 
coupling of the time-dependent metric component g+ + to the massless field S. 

In order to compare the matrix model with the spacetime results we must know 
whether the matrix model result corresponds to conformal, inertial, or yet another coordi- 
nate choice. Although the literature is not very explicit on this point, in most discussions of 
the matrix model it is implicitly assumed that the matrix model time t mrn and eigenvalue 
are related to the spacetime retarded time via 

t mm - In \ = tr nerUal . (7.4) 

This assumption has been adopted in this paper in equation ( |3.1|) , and we have seen herein 
that the full string theory answer is consistent with the matrix model answer only when 
we make this (standard) assumption. 

In this paper we are considering normalizable excitations of the ground state. In the 
context of large, non-normalizable perturbations, it is less obvious how time should be 
defined. Some discussion is in [|il|[T4] , |l6|| . 
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7.2. The non-singlet sector 

In this paper we have analyzed whether black holes could form in the singlet sector. 
Critics have pointed out that our efforts were doomed to fail because it is known that 
the SL(2)/U(1) black hole involves non singlets. However, it is also known that the 
SL(2)/U(1) solution is not the only one that looks like a black hokH. So it might well be 
that there are other solutions, which exist purely in the singlet sector, which also have a 
black hole interpretation, even if they are not the standard SL(2)/U(1) black hole. Our 
results, however, are consistent with the idea that there are no black holes in the singlet 
sector. So a plausible explanation for the peculiar horizon-aversion of collapsing shells 
found in section 3 is that a black hole cannot be formed because the matter is in the 
singlet sector and the black hole is not. 

If the non-singlet nature of black holes is the explanation, one may also ask if black 
holes can be formed from non-singlet initial data. For example one might try to form them 
by collapsing DO-branes, which also can be argued to involve non-singlets. It would also 
be nice to know what the matrix model dual to a Lorentzian black hole is, as opposed to 
euclidean black hole ||. 

7.3. SL(2) k /U(l)=black hole for small k? 

For large k, Witten's SL{2)k/U{l) coset is a weakly coupled sigma model describing 
a black hole. For small k the sigma model becomes strongly coupled so it is no longer clear 
if it should be thought of as a black hole. In the matrix model context, the relevant cosets 
always are in the small k regime. Therefore it is not clear if they should be interpreted as 
black holes, as we discuss in this section. 

As it was discussed in |^9|] the nature of the SL{2)k/U{l) confer mal field theory de- 



pends in an important way on k. For large k, k > 3, this confer mal field theory has a 
normalizable zero mode which corresponds to changing the value of the dilaton at the 
tip0. For k < 3 this mode becomes non-normalizable. This can be seen from the alge- 
braic analysis of the SL(2)k/U(l) model as follows. The zero mode is constructed from 
discrete representations as JZ\ J~\ \j, Tn) with j = m = lEl. This is a state with worldsheet 



16 There are configurations that look like black holes with different radii []5| . 

17 Actually, the theory has two zero modes one changes the value of the dilaton and the other 
changes the integral of Bns on the cigar. 

18 In conventions as in [ 30 1 . See [31| for further discussions about this state. 
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conformal weight one. In general, a state is normalizable only if ~ < j < ^ — |, so this 
state with j = 1 is normalizable only for /c > 3. From the target space point of view this 
seems surprising. It is clear that the wavefunction of the vertex operator that changes 
the dilaton at the tip decays as e~ 2 ^ for <fi — ► oo while the normalizability condition only 
requires that it decays faster than 0. The crucial point is that a mode on the cigar 
does not have definite winding at infinity. A good way to think about it is to imagine that 
the wavefunction of a given mode contains all possible winding numbers. In particular, 
the wavefunction of this mode goes as 

V s * ~ JZ x JZ x e~^ + e -2(fc/2-i)* e Wfc(x*-*«) + e -2(k/2-i)d, e -iSk(x L -x R ) + . . . (7 5) 

where the dots indicate higher winding modes. We see that the winding number one terms 
in ( [7.5| ) are not normalizable for k < 3. 

The conclusion from these remarks is that the cigar conformal field theory at k < 3 
does not have a zero mode and furthermore, when we go to the weakly coupled region, it 
differs from the linear dilaton theory by a non-normalizable operator. This operator is the 
one appearing in the second term in ( [7.5| ) . Note that we have derived all this from analysis 
purely in the SL(2, R)k/U(l) theory. It is not necessary to assume that there is a dual 
sine-Liouville theory. Of course, the fact that the operator is not-normalizable is clear in 
the sine-Liouville theory and in a different context it was noted in ||32|1 . Even though this 



can be derived purely from the SL(2, R)k/U(l) CFT, these facts cannot be derived from 
the lowest order in a' analysis of the sigma model, even if we use the supposedly "exact" 
metrics in the literature. Any analysis that focuses only on the lowest order in a' will miss 
the winding modes, which are not included. The point is that the curvatures are large and 
it is necessary to solve the theory exactly in a' . 

For k > 3 the existence of this zero mode is conceptually related to the fact that the 
asymptotic value of the radius of the cigar cannot be changed, any attempt to change the 
value of the radius will drive the value of the dilaton at the tip to zero or infinity. On the 
other hand, for k < 3 there is no zero mode, so we expect that it should be possible to 
change the value of the radius. Indeed, the matrix model analysis of [|5) shows that the 
value of the radius can be changed. 
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The last statement is saying how we are normalizing cf>. 
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All of these points are highlighting the fact that the Lorentzian continuation of the 
cigar conformal field theory is not obvious. In particular, we would need to understand the 
Lorentzian interpretation of the extra winding mode operator that is being turned on@. 

One can, nevertheless, be very naive and compute the absorption cross section for the 
black hole. If the absorption is nonzero it is reasonable to call the geometry a black hole. 
For this one can compute the momentum two point function in the cigar geometry and 
analytically continue the result to Lorentzian signature. This gives a non-zero answer for 



the absorption p^][pl[, even when k < 3. If we believe this continuation it can be taken 



as evidence that the SL(2, R)k/U(l) CFT is a black hole even for k < 3. 
7.4- Conclusions 

Previous investigations of high energy tachyon scattering 00], have shown that it 
does not lead to the formation of black holes. In those studies scattering of a single high 
energy tachyon was considered. In this paper we considered the collapse of a coherent 
state formed by many tachyons. By looking at the time at which the energy is emitted 
we concluded that there cannot be a long lived black hole. The energy is emitted within 
a time of order log E of the the time at which the black hole would be formed. If a black 
hole had formed we would have expected that a finite fraction of the incoming energy 
would come out over a time proportional to E. This simple argument does not rule out 
the formation of a black hole with very high temperature. 

On a different, but related note, we have also asked here the following question: to 
what extent is the low energy effective action corrected? We know from the analysis of |23[] , 



that many of its features are indeed consistent with the matrix model. In particular, [p3 | 
have shown that the matrix model includes gravitational interactions which are present 
in the linear dilaton theory, in the region far from the Liouville wall. Particle creation is 
another effect mediated by the gravitational field, so our naive expectation was that we 
were going to see it if we repeated the analysis of [E3|. However, here there is an interesting 
cancellation between the gravitational effect and the effects coming from the exchange of 
all massive fields. This cancellation, which is present in the exact string theory answer, 
implies that there is no net particle creation by an infalling pulse. 



20 A naive Lorentzian's continuation of the matrix model in [|| is not possible because it would 
lead to a Lagrangian that is not local in time, since the Wilson lines are in the exponent. Maybe 
a local Hamiltonian might be obtained after "integrating in" some quark degrees of freedom... 
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Our computation of particle creation by the infalling pulse focused just on the leading 
term. It would be nice to see what happens at higher orders in the e* expansion. 
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